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theoretical value of the spectral index and its running compatible with the latest PLANCK data, 
to have a reheating process via gravitational particle production, or to predict some signatures in 
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1. Introduction 

It is well-known that inflation suffers from several problems (reviewed for instance in [|I|]), like 
the initial singularity which is usually ignored, or the fine-tuning of the degree of flatness required 
for the potential in order to achieve successful inflation [^. 

An alternative scenario to the inflationary paradigm, called Matter Bounce Scenario (MBS), 
has been developed in order to explain the evolution of our Universe |^] avoiding those problems. 
Essentially, it presents at very early times a matter dominated Universe in a contracting phase, 
evolving towards the bounce after which it enters an expanding phase. This model, like inflation, 
solves the horizon problem that appears in General Relativity (GR) and improves the flatness prob¬ 
lem in GR (where spatial flatness is an unstable fixed poinf and fine funing of inifial condifions is 
required), because fhe confribufion of fhe spafial curvafure decreases in fhe confracfing phase af fhe 
same rafe as if increases in fhe expanding one (see for insfance |0]). 

The aim of our work is fo show viable bouncing cosmologies where fhe maffer parf of fhe La- 
grangian is composed of a single scalar held and, fherefore, have fo go beyond General Relafivify, 
since GR for flaf spafial geomefries forbids bounces when one deals wifh a single held. Hence, 
fheories such as holonomy correcfed Loop Quanfum Cosmology (LQC) |^, where a big bounce 
appears owing fo fhe discrefe sfrucfure of space-fime |^] or feleparalellism, fhaf is F{T) gravify, 0 
musf be faken info accounf. 

The unifs used in fhe paper are: h = c = SnG = 1. 

2. F{T) gravity in flat FLRW geometry 

Teleparallel fheories are based in fhe Weitzenbock space-time. This space is wifh a Lorenfz 
mefric, in which a global, orfhonormal basis of ifs fangenf bundle given by four vecfor fields {e,} 
has been selecfed, fhaf is, fhey salisfy g{ei,ej) = rjij wifh rj = diag (—1,1,1,!). The Weifzenbock 
connecfion V is defined by imposing fhaf fhe basis vectors e, be absolufely parallel, i.e. fhaf Ve, = 0. 

The Weifzenbdck connecfion is compafible wifh fhe mefric g, and if has zero curvafure be¬ 
cause of fhe global parallel fransporf defined by fhe basis {e,}. The information of fhe Weifzen¬ 
bock connecfion is carried by ifs forsion, and ifs basic invarianf is fhe scalar torsion T. The 
connecfion, and ifs forsion, depend on fhe choice of orfhonormal basis {e,}, buf if one adopfs 
fhe flaf Friedmann-Lemaifre-Roberfson-Walker (FLRW) mefric and selecfs as orfhonormal basis 
{eo = do,e\ = ^di,e 2 = ^^ 2,03 = ^^ 3 }, then fhe scalar forsion is T = — 6 H^, where H = ^ is fhe 
Hubble paramefer, and fhis identify is invarianf wifh respecf fo local Lorenfz fransformafions fhaf 
only depend on fhe time, i.e. of fhe form e, = K\{t)ek (see [^, ^). 

Wifh fhe above choice of orfhonormal fields, fhe Lagrangian of fhe F{T) fheory of gravify is 

^T = r(F(r)+j2’M), (2.1) 

where F = a^ is fhe elemenf of volume, and .^m is the matter Lagrangian density. 

The Hamiltonian of the system is 

M’T={lT^^^-F{T)+p^'r, (2.2) 
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where p is the energy density. Imposing the Hamiltonian constrain = 0 leads to the modified 
Friedmann equation 


p = -2^^T + F{T) = G{T) 


(2.3) 


which, as r = —6H^, defines a curve in fhe plane {H,p). 

Equation (^) may be inverfed, so a curve of fhe form p = G{T) defines wF{T) fheory wifh 


F{T) = - 


f G{T) 


■/ 


T^/^ 


dT. 


(2.4) 


To produce a cyclically evolving Universe, lef us fake fhe F{T) fheory arising from fhe ellipse 
fhaf defines fhe holonomy correcfed Friedmann equafion in Loop Quanfum Cosmology 




(2.5) 


where Pc is fhe so-called critical density. 

To obfain a paramefrizafion of fhe form p = G{T), fhe curve has fo be splif in fwo branches 


p = G±(r) = ^(i±Ji + ^ 


( 2 . 6 ) 


where fhe branch p = G-{T) corresponds fo // < 0 and p = G+(r) is fhe branch wifh H > 0. 
Applying Eq. (2.4) fo fhese branches produces fhe model (|^ [I^]) 


Tpc 


F±{T) = ±\l —^arcsin 


2T 

Pc 


+ G±(r). 


(2.7) 


3. Matter Bounce Scenario 

Mailer Bounce Scenarios (MBS) [||] are essentially characferized by fhe Universe being nearly 
mailer dominated al very early limes in fhe conlracling phase (fo obfain an approximalely scale 
invarianl power spectrum), and evolving towards a bounce where all the parts of the Universe 
become in causal contact solving the horizon problem, to enter into a expanding regime, 
where it matches the behavior of the standard hot Friedmann Universe. They constitute a promising 
alternative to the inflationary paradigm. 

According to the current observational data, in order to obtain a viable MBS model, the bounc¬ 
ing model has to satisfy some conditions that we have summarized as follows: 

1. The latest Planck data constrain the value of the spectral index for scalar perturbations and 
its running, namely and to 0.9603 ± 0.0073 and —0.0134 ± 0.009 respectively []ri|]. 
The analysis of these parameters provided by Planck makes no slow roll approximation (for 
example, the first year WMAP observations was done considering the ACDM model [|T2[]), 
and thus, the parameters and can be used to test bouncing models. A nearly scale 
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invariant power spectrum of perturbations with running is obtained either as a quasi de Sitter 
phase in the expanding phase or as a nearly matter domination phase at early times, in the 
contracting phase [jl^. Then, since for an exact matter domination in the MBS the power 
spectrum is flat, to improve the model so as to match correctly with the observational data, 
one has to consider, at early times in the contracting phase, a quasi-matter domination period 


characterized by the condition 
and the energy density of the Universe. 




<C 1, with P and p being respectively the pressure 


2. The Universe has to reheat creating light particles that will thermalize matching with a hot 
Friedmann Universe. Reheating could be produced due to the gravitational particle creation 


in an expanding Universe []14|]. In this case, an abrupt phase transition (a non adiabatic 
transition) is needed in order to obtain sufficient particle creation that thermalizes producing 
a reheating temperature that fits well with current observations. It is shown in [ ]I5| ] that 
gravitational particle production could be applied to the MBS, assuming a phase transition 
from the matter domination to an ekpyrotic phase, which also maintains the isotropy of 
the bounce, i.e., solves the Belinsky-Khalatnikov-Lifshitz instability [|T^, in the contracting 
regime. In [[T^, dealing with massless nearly conformally coupled particles, a reheating 
temperature compatible with current data has recently been obtained. 

3. The data of the seven-year survey WMAP ([|T8|]) constrains the value of the power spectrum 
for scalar perturbations to be ${!<■) — 2 x 10^^. The theoretical results calculated with 
bouncing models have to match with that observational data. 

4. The constrain of the tensor/scalar ratio provided by WMAP and Planck projects (r < 0.11) is 
obtained indirectly assuming the consistency slow roll relation r = 16£ (where e = —jjj = 

A (^') is the main slow roll parameter) [p^, because gravitational waves turn out not to 


have been detected by those projects. This means that the slow roll inflationary models must 
satisfy this constrain, but not the bouncing ones, where there is not any consistency relation. 
Actually, to check if the MBS models provide a viable value of the tensor/scalar ratio, first 
of all gravitational waves must be clearly detected in order to determine the observed value 
of this ratio. 


5. Some non-gaussianities has been detected by PLANCK’s team [|20|]. A theoretical viable 
bouncing model has to take into account these non-gaussianities. However, as shown by the 
recent calculation of the 3-point function in bouncing cosmologies this seems to be a 
very difficult challenge in MBS. 


4. Perturbations in Matter Bounce Scenario 


The Mukhanov-Sasaki equations for F(r) gravity and LQC are given by [p2|, p3|] 



zi 


-cy^i;siT) + ^QT) = o 


'S{T) 


(4.1) 
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where and denote the amplitude for scalar and tensor perturbations. 
In teleparallel F{T) gravity one has [ ^ ] 


Z<: = 




■\2 


cW 


Zt = 


2 2 
a ct 


arcsin (2 


inl 




H 


with = 1 — 


2p 


(4.2) 


H 


In contrast, for holonomy corrected LQC one has [25] 

-,2 


Zs = 


2 • 2 


Zt = 


a 


ci = a. 


(4.3) 


a’ 

Dealing with the simplest model of MBS, i.e., a background depicted by an exact matter 
domination, whose scale factor is given by a{t) = (|pct^ + l) 
spectrum for scalar perturbations, given by 

2 


one obtains the following power 


3p^ 

^s{k) = 

Ppi 


Zs\ri)dri 


(4.4) 


In the particular case of teleparalell F{T) gravity one has =^s(k) = 


Q where ^ = 

I — ^ + jj — ... = 0.915965... is Catalan’s constant, and for holonomy corrected LQC 3^s{k) = 


E-Bl. ||2^ 

9 Pp, 

In the same way, the ratio of tensor to scalar perturbations in MBS is given by 


r = 


8 r^z-j 




3 \ f“„Z/(r])dr] j 


(4.5) 


obtaining for teleparallel F{T) gravity r = 24 ) ’ '''here Si{x) = /q ^^dy is the Sine inte¬ 

gral function. In contrast, for holonomy corrected LQC this ratio is very close to zero. 

A final remark is in order: All of the current models that depict the MBS have some problem. 
For example, in holonomy corrected LQC the square of the velocity of sound becomes negative in 
the super-inflationary phase (Pc/^ < P ^ Pc)- spite of the fact that in holonomy corrected LQC, 
in order to calculate the power spectrum, only modes that satisfy the long-wavelength condition 
are used, it is important to realize that, in the super-inflationary regime, the ultra- 


c?k2| <C 


violet modes (modes that satisfy the condition \clk^\ S> 


) will suffer Jeans instabilities, and 


thus undesirable cosmological consequences could appear. This is a serious problem that needs 
to be addressed in holonomy corrected LQC. In F{T) gravity, the square of the velocity of sound 
is always positive. However, teleparallelism suffers from the problem that the main invariant, 
the scalar torsion T, is not at all an invariant, in the sense that it depends on the choice of the 
orthonormal basis in the tangent bundle (the tetrad). The scalar curvature R is a real invariant, 
but the current bouncing models in modified F(R) gravity, do not support a matter domination in 
the contracting phase, and thus their power spectrum is not nearly scale invariant /|2^. 

Another possibility is to consider the MBS for a flat FLRW geometry in the context ofGR. In 
this case one needs more than one field, one of them a ghost condensate field ora Galilean type 

field /|5^ which violates the Null Energy Condition, to obtain a non-singular bounce. In this case 
the scenario is very complicated, and it is not clear at all how to compute the theoretical quantities 
such as the spectral index and its running, to be compared with the current observational data. 
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